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The flexible nature of exchange, correlation and Hartree physics: resolving 
"delocalization" errors in a 'correlation free' density functional 

Tim Gould-"^ and John F. Dobson-"^ 

Qld Micro- and Nanotechnology Centre, Grijjith University, Nathan, Qld 4^1^! 
Australia 

By exploiting freedoms in the definitions of 'correlation', 'exchange' and 'Hartree' physics in ensemble systems 
we better generalise the notion of 'exact exchange' (EXX) to systems with fractional occupations functions 
of the frontier orbitals, arising in the dissociation limit of some molecules. We introduce the Linear EXX 
("LEXX") theory whose pair distribution and energy are explicitly piecewise linear in the occupations f,^. 
We provide explicit expressions for these functions for frontier s and p shells. Used in an optimised effective 
potential (OEP) approach it yields energies bounded by the piecewise linear 'ensemble EXX' (EEXX) energy 
and standard fractional optimised EXX energy: iJ^EXX ^ ^lexx ^ £'EXX_ Analysis of the LEXX explains 
the success of standard OEP methods for diatoms at large spacing, and why they can fail when both spins 
are allowed to be non-integer so that "ghost" Hartree interactions appear between opposite spin electrons in 
the usual formula. The energy ijLEXX contains a cancellation term for the spin ghost case. It is evaluated for 
H, Li and Na fractional ions with clear derivative discontinuities for all cases. The p-shell form reproduces 
accurate correlation-free energies of B-F and Al-Cl. We further test LEXX plus correlation energy calculations 
on fractional ions of C and F and again shows both derivative discontinuities and good agreement with exact 
results. 

PACS numbers: 31.15.ep,31.15.eg,31.10.-fz 



I. INTRODUCTION 

Following initial work by Yang and coauthors^""* on 
non-interacting ensembles^ with spin-resolved fractional 
occupancy, much consideration has been given to the be- 
haviour of density functional theory (DFT) under the 
Kohn-Sham (KS) prescription^, and its various common 
approximations (eg. LDA*^, GGA^, Becke-like^, OEP^) 
in such ensembles. Many attempts have been made to 
understand and deal with the issues that arise in ensem- 
bles (see eg. Refs 10-12), with variable success. We 
will show that, in such systems, the notion of 'correla- 
tion' physics becomes intertwined with 'exchange' and 
'Hartree' physics in the usual prescription, with (improv- 
able) consequences for common approximations. 

Let us begin by considering, quite generally, the na- 
ture of 'electron correlation' and 'electron exchange' in 
a non-ensemble system. The usual expression for the 
groundstate correlation energy can be written as 

E^^ i'^lHl^) ~ {^^\H\^'^) (1) 

where H is the Hamiltonian of a many-electron system, 
4*) is its groundstate wavefunction, and |^'"^) is some 
approximation to the wavefunction (by the variational 
principle, correlation energy is never positive). Thus cor- 
relation is not an intrinsic property of the system, but a 
property of the chosen trial wavefunction. In standard 
optimised effective potential (OEP) approaches^, includ- 
ing KS DFT, 1^"^) takes the form of a single Hartree- 
Fock like Slater determinant which is constructed from 
one-particle orbitals \i(7) evaluated in a common one- 
particle Hamiltonian h = i + V^^ where t = — ^V^ and 
V = Vcr{r). We can now define the exchange energy = 
(^'^|//|\E''^) — E and the "naive Hartree" energy of the 



system" E = + '^ExtN^T) + U ^n{r)nir'). 

Here n{r) = (vf"^ |^('")| ^''^) = J2ia- l'/'io'('')P [where h{r) 
is the electron number density operator and (piair) = 
{r\ia)] and Vexi = ^Ext('") is the external potential. The 
groundstate energy is thus E = E + E^ + E^ where the 
partitioning depends on both the choice of E and l^*^}. 

This can be extended into ensembles by replacing pro- 
jections on wavefunctions O = (4'|0|^) by traces on den- 
sity matrices O = Tr[/50] (where operators act appropri- 
ately for any number of electrons) and by summing E 
over ensemble members. The density matrix p is defined 
as 



(2) 



where < W£ < 1 is the weight of member £ with wave- 
function and ^£W£ — 1. Minimisations can then 
be carried out over p rather than |<I>). 



II. EXACT EXCHANGE APPROACHES 

We can now succinctly define the standard 'exact ex- 
change' (EXX) functional approach. Here we consider 
only E^^-^ — E + Ex with Ec assumed to be zero. In- 
vestigations into EXX in fractionally occupied ensem- 
ble systems^'^^"^^ show both successes and shortcomings 
(discussed in more detail later). In all these works, the 
Hartree and exchange energy takes the 'standard' form, 
bilinear in the occupations f[: 



drdr' 
2|r - r' 
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where Piaja' = \(l)ta{r)\'^\(t)]a'{r')'^\ and Qiaja = 
(^icr(^)0*cr(^')'/'^cr(^)'/'jo'(^')- Hcrc thc negative exchange 
term cancels thc unphysical positive Hartree interaction 
of each spin orbital \ia) with itself. However if two 
different orbitals of the same spin are partly occupied 
(0 < fi,fj < 1 with i ^ j), or if there is partial occupa- 
tion of both spins in thc same orbital (0 < //,/^ < 1), 
there is a corresponding cross-term in (3) that is not can- 
celled. 

In a slightly different context Gidopoulos et al}^ call 
this spurious term the "ghost interaction" as it represents 
an unphysical interaction bcrt.wccn orbitals in different 
non-interacting ensemble members. In the regular EXX 
energy expression (3), the ghost interaction appears in 
the Hartree and exchange energy terms involving pairs 
of orbitals in the frontier orbital. In a Kohn-Sham in- 
terpretation of the equivalent diatom problem, these in- 
teractions would be supressed in the total energy via or- 
thogonality of the degenerate groundstate wavefunctions. 
However, when one does not have the exact exchange- 
correlation functional, or as here neglects correlation, it 
can reappear, particularly when one does not properly 
account for the ensemble nature of the system. 

We will argue that, in the ensemble interpretation 
of partial occupation^"^, this cross term should not be 
present, and its explicit removal results in an improved 
linear exact exchange (LEXX) approach which is cor- 
rectly piecewise linear, not bilinear, in the occupation 
factors /. Here, defining to be one for orbital jzcr) 
occupied in ensemble member £ and zero otherwise, 
we exploit the fact that the 'ensemble occupancy' fac- 
tor f^ = {ODs = Ylic requires weights ws that 
are piecewise linear in /f , from which it follows that 
{9'^6'j )s = J2e '^£^is^js similarly piecewise linear. All 
energy terms arc proportional to {0f)£ or {BfOJ and 
are thus piecewise linear. As will be discussed in more de- 
tail later this is equivalent, under an exchange approach, 
to finding a non-interacting ensemble of Slater determi- 
nants formed from a common set of orbitals produced in 
a common potential. 

This allows the creation of simple functional that 
avoid much of the "localization and delocalization error" 
of Yang ct al.^^*, and thc "many electron self interaction 
error" of Perdew et al.^®. In the present work we focus 
on two illustrative cases: i) a single partially occupied 
"frontier" orbital with < < 1 and < < 1; 
and ii) open p shells with fj^=fj^. However the scheme 
itself has wider applicability, including the full dissocia- 
tion problem of molecules. Ref. 18 might be considered 
another specific example of this approach, while Ref. 20 
outlines a similar approach via HF for the restricted case 
of fractional occupation of a single spin (their ISSO ap- 
proach) . 



A. Non-interacting 'exchange' ensembles 

To illustrate the general approach we consider, as an 

example, ensembles with total and spin-resolved electron 
number Nt = N + f and Nt„ = N/2 + f {N is even). 
The groundstate ensemble members and weights can be 
found be minimising over density matrices subject to var- 
ious constraints. However for simple cases where energy 
ordering is obvious, one can construct the ensemble more 
intuitively, just by demanding that a given set of occu- 
pations ff be reproduced. For example if the frontier 
orbital is non-degenerate (eg. in an s shell), then the en- 
semble will be composed of up to three components. For 
/ < 1, the ensemble is formed from parts an A'^-l-l elec- 
tron system with extra electron in ^ (short-hand N+ t) i 
parts A^-|- 4- and (1 — /) parts where, because A^ is 
even, both spins are filled equally. For / > 1 the ensem- 
ble comprises (1 — /■'') parts A'-l- t, (1 — /''") parts N+ 4-, 
and (/ - 1) parts jV + 2. 

The density matrix is composed of many-electron 
wavefunctions |3>f ) and is 

= Y^W£\^£){<^e\- (4) 

s 

For the present case of a non-degenerate frontier orbial 
W£ & {I - f,f,f^} and $f e {$jv,$Ar+-t-,$Ar+4} for 
/ < 1 while W£ & {I- /^l - f\f ^ 1} and S 
{<&Ar+-t-, $Ar+4., $Ar+2} for / > 1. This leads to a total 
energy E{f) = Tr[/3-^ H] = Y,s W£E[^s] that obeys 

rr(f.UEM+l + {l-f)EM, 0</<l 

\{f-l)EM+2 + {2-f)EM+i, l<f<2^' 

where Ef^ is the energy of an A-electron system (note 
that i?7v+t — ^N+i = -Eat+i). 

The LEXX is defined in general by assuming that the 
trial density matrix p^^ of the ensemble obeys the same 
relationship (4) h\it with the component wavefunctions 
l^f) replaced by Hartree-Fock like determinants |^>£) 
constructed from a single set of spin-dependent orbitals 
{jio")}. This trial density matrix: i) reduces to the regu- 
lar EXX for integer occupation, ii) gives correct energies 
for H with less than one electron, split arbitarily between 
spins, and iii) is constructed from a single set of orbitals 
\ia) evaluated in a common Hamiltonian, a requirement 
that ensures that OEP or KS methods can be used. Here 
the orbitals are eigen-solutions h\i(j) = eia\ia) of a one- 
body Hamiltonian h = t + V . We sort the orbitals so 
that EiCT < Ejcr for i < j. Taking the spin-resolved density 
na{r) = Tv\p^'^ ha{r)] one now finds 

n.{r) =Y,{ei)e\cl>ia{r)\^ ^ /f |0,.(r)|^ (6) 

i i 

where typically f^ = 1 for the inner orbitals and = 
where \ha) is the frontier orbital in the spin-shell with 
highest energy: which may or may not be occupied in 
both spins. 
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The EXX approximation {E^ = 0) allows us to use 
only the Hartrce and exchange (Hx) components of the 
pair-density n2Hxcrcr' = T^r[pf'^ha{'r')ha'{r')] to evaluate 
the electronic groundstate. From the properties of HF 
wavefunctions, the pair-density of an ensemble can be 
written as 



= ^(^f^J )e[Pi<Tja' — Saa'Qitrjtr]- 

ij 

Finally, we can use (7) to calculate the energy via 
1 /' drdr' , 

+ 2^ J h;r-;^"2Hx.a'(r,r) 

(7(7' 

where t^r) = ^|V0i.(r)|2 and 



(7) 



(8) 
(9) 



el 



(1) _ 



= / dr 



icrjo 



1 f drdr' 

cr' ~2 J 1^ _ ^/| [^icr3(^' ~ ^cra'Qiaja] 



crcr'^i<Tj<T\ • 



(10) 

(11) 



These energy expression are perhaps the most general, 
and most important in this work, highlighting the impor- 
tance of ensemble averages in the evaluation of average 
occupation and pair-occupation factors for groundstate 
energy calculations. 



B. Fractional s shells 

For the fractionally occupied s shells discussed here, 
{e^0f}s = min[/f , //'] - S.hjhSaa'C^ (a is the opposite 
spin to a and Cj} is defined below). The Hartree and 
exchange components can be compactly written as 

n2Haa' = J2 "^Hfi, ff]Pi^j^' ' ^aa'C^Phaha, (12) 

ij 

n2^aa' = - Saa' ^ ^Hfi . fj]Qiaja (13) 



where we have chosen to split Hartree and exchange 
terms via P and Q. The term 



C^=mm[f\f\{l-f),{l-f^)] 



(14) 



removes spurious "ghost interactions" between electrons 
of unlike spin. For a zero to two electron system equa- 
tions (12)-(14) are equivalent (after integration) to equa- 
tion 7 of Ref. 4 sans the correlation energy term. This 
desirable outcome is a direct result of the ensemble aver- 
aging. 



When either or /-•- is integer, C^} = and n2Hxaa' = 

^^ij fi fj [Piaja' ^aa' Qiajcr'] since Piaia Qiaitj' 

Clearly this is the form used in (3) and thus energies de- 
rived from (12)-(13) will be identical. We can now proffer 
an explanation for the variable success of the EXX for 
fractionally occupied ensembles. By violating the aufbau 
principle and/or allowing spins to vary in an unrestricted 
fashion, good results have been obtained for atoms and 
diatoms^^'^^ and systems with fractional occupancy^. In 
these works only one spin was allowed to be non-integer 
so that {0'^e'^')£ = f^fi and the EXX and LEXX ener- 
gies were equivalent. In systems where both spins were 
fractionally occupied (eg. Refs. 15 and 17) the EXX 
failed to reproduce the correct derivative discontinuity. 
In these works = = //2 and ((?^(?^')£ + fj^f;:' . 
Thus the EXX and LEXX energies differed. We show 
later that, in this case, the LEXX is guaranteed to pro- 
duce a lower energy. 



C. Fractional p shells 

As a less trivial example, we also consider the case 
of degenerate frontier p orbitals with equal densities in 
each spin. Here we must sum not only over ensembles 
members of different electron number, but also over the 
degenerate combinations of p^, py and pz orbitals. For 
example, in an isolated carbon atom each ensemble mem- 
ber has fully occupied Is and 2s shells, but only two 
occupied 2p orbitals of the same spin a which we de- 
note p^a and psa where j ^ 6 and 7,(5 G {x,y,z}. To 
find the equal-spin, spherically symmetric ensemble we 
weight each ensemble equally so that Wp^p^a = | for all 
six combinations of 7 ^ 5 and a. In member p-yPsa we 
set = ^2p P5 ~ 1 while the remaining 2p orbital 

with spin cr, and all 2p orbitals with spin tr have zero oc- 
cupation. Averaging over all cases gives {02p,p^)£ — ^ ^ 
expected, while {0^p^p_0^p^pj£ = |, {02p,p^0^p,ps)£ = g 



for 7 ^ (5 and 



= 0. 



For general unfilled frontier p shells this yields an ad- 
ditional like-spin correction of the form —C^"'[P — Q] to 
(12) and (13) so that 

ij 

- {5aa'Cl - 5„^,C^) Phaha', (15) 
h 

n2xaa' = - E min[/f , fj]Qiaja 

ij 

-S^^^C^Y^Qh^ha, (16) 



where we recognise the degeneracy in the outermost p 
orbitals by summing over h with equal weights. Let us 

restrict ourselves to the case = /■'' = //2 where < 
/ < 2 is the total occupation (over both spins) of each 
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orbital in the shell. One can sum over the ensemble to 
show (after much work) 

C^=imin[/,2-/,|l-/|,l/3] (17) 

for open p shells. We note that the total number of elec- 
trons in the shell is Np = 3/ and (15)-(17) are valid for 
Np integer or fractional. 

The like-spin correction ensures that a bilinear ap- 
proach would fail even for systems with one spin fully oc- 
cupied. Indeed it is only true for the case = 1, = 
(or vice versa) occurring for half-occupied shells in N and 
P. Here one must not only allow the spin-symmetry to be 
broken, but also break the spherical symmetry to make 
the bilinear expression (3) correct. 

D. General ensemble systems 

While we have so far determined our ensembles using 

explicit knowledge of the degenerate groundstate, it is 
possible to carry out a more general ensemble minimi- 
sation to determine wg. Here, for a given potential, we 
allow ensemble members (determined by member occu- 
pancy factors ^^) to sample all combinations of 'occu- 
pied' and 'unoccupied' orbitals of the one-electron Hamil- 
tonian, and minimise the energy with respect to ws- In 
practice we would restrict the allowed ensemble members 
to limited combinations predicted to be low in energy. 
For example in the p shell case given above, or indeed 
Be, we might search for the minimum over cases with 
full occupancy in Is^ and varying occupancy in the near- 
degenerate 2p and 2s orbitals. 

As shown in equation (9), the general LEXX energy 
^LEXxji^^jj £qj. ^ given ensemble can be written as 
a sum of the ensemble averaged occupations {Of)£ and 
pair-occupations {6f6j )s with orbital dependent energy 
prefactors given in equations (10) and (11). These aver- 
aged occupations depend piecewise linearly on ws and 
thus can be minimised under the con- 

straints < < 1 and ^£ W£ = 1 ie. we look for 
the (constrained) set of weights minimising 

^LExx =^^,[^^^,e« + O^eeie^,^- (18) 

£ ia iaja' 

We can thus find, for a given potential and orbtials, 

the optimal weights wg. and through them (^f)^ and 
(^r^J )f • For the true KS potential, this should be equiv- 
alent to finding the temperature— > 0+ limit of finite- 
temperature DFT. Such an approach might be useful for 
dealing with the difficult atomic dissociation problem. 

III. OPTIMISED EFFECTIVE POTENTIALS 

For a many-electron system the EXX (or LEXX) 
groundstate energy is composed of the orbital kinetic en- 



ergy = i / dr Yjia ./r|V(/)icrP, the energy from the ex- 
ternal potential -Esxt = / drVExt^ and the Hartree plus 
exchange energy £^hx- For an ensemble we calculate E'hx 
via, for example, the expansion (12)-(13) of n^Hi^aa' for 
s shells [or (15)-(16) for equi-p shells] to form the orbital 
dependent LEXX expression 

E^-=Y.I ^?^«2Hx..'(r,r'), (19) 

while for 'standard' EXX we instead use (3). The differ- 
ence in energies between the LEXX and 'standard' EXX 
for frontier s shells is thus the difference between (19) 
and (3). E.g. for the s case 

^LEXX _ ^EXX ^ _ ^ _^h^^ (20) 

where is given by (3) and en = J ^^Phm and 

Cu=Cu- mini/l', /^] + f\f^ = min[/t/i, (1 - /t)(l - 
/"'')] governs the unlike-spin correction to the Hartree en- 
ergy required when both f"^ and /■'- are non-integer. A 
similar expression can be derived for the like-spin correc- 
tion to p-shells. 

We can now define orbital dependent groundstate en- 
ergies via E^^^ = T,+ Ee^i + E^^ for the EXX and 
^LEXX ^ ^Exx _ (jh^^ fQj. ^j^g LEXX. In an optimised- 
effective potential^ approach, we look for a potential 
V = Vocr{r) such that the orbitals satisfying [i+Voa]4'icr = 
eia4>icr minimise the energy. Here we call this approach 
the EXX or LEXX (with an overline to denote use of an 
optimised effective potential) depending on the Hx func- 
tional used. Finding V^^ involves, as input, the functional 
derivatives Dia-{r) = SE-ii^/S(j)ia-(r). Thus the scheme for 
finding optimised LEXX solutions differs only from that 
for the regular EXX in that Djo- for the LEXX includes 
an extra term for i = h. Via C^j, the additional term 
vanishes whenever f"^ or /■'' is integer, as expected (at 
least for s shells). 

Let us consider some of the formal implications of the 
LEXX. Firstly, the total energy found in an optimised 
LEXX scheme must be bounded below by the EXX en- 
ergy of the full ensemble. To prove this we first note 
that the ensemble EXX energy £;eexx f^j. g^j^ g^semble 
of positive weights W£ of elements £ can be written as 
£;™(/) = Eey^sEf^^MJ] where [t + ^lC = 
ef^(/>f^ and V^^ is chosen to minimise and may 

vary between different ensemble members. From (6)- 
(19), it is clear that i;LExx||^.^|] = j^e «;f£;f^^ 
where Voa in [t + Voct]4> ia- — ^ia4'ia Can no longer vary 
separately for each part of the ensemble. Thus by the 
variational nature of an OEP we find E^^^[{4>^^}] < 
Ef'^'^licpia}] and i;EEXX(y) < ^lexx Secondly, we see 
that E^^^^[{(j)}] < E^^^[{(f>}] for any set of orbitals {<!>} 
and thus £;LEXX[LEXX] < E^^^^[EXX\ < E^^^[EXX] 
(where the term in the square brackets labels the OEP 
used to evaluate the orbitals) with the equality hold- 
ing (for s shells) only when Cjj = (ie. when each 
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of the spins is integer occupied). The former inequal- 
ity follows from (20) by noting that > and e'* = 
/ ^^\Ph-thi > as Ph^hi > (similarly for the like spin 
term) and the latter follows from the minimisation prin- 
ciple of OEPs. Putting the OEP inequalities together, 
we find 



FIG. 



Groundstate energy differences E{f,f^) 



(21) 



where we include the overline (indicating an optimised 
potential was used) for clarity. 



IV. CORRELATION ENERGIES 

The consequences of the improved pair-densities also 
extends beyond exchange physics. Some beyond-dRPA 
correlation energy methods [see Ref. 21 for an overview] 
like the RPAx22, RXH^^ and PGG^^ kernels, ISTLS^s 
and tdEXX^^ depend in some way on the groundstate 
pair-density. The difference between the EXX and LEXX 
expressions will therefore manifest in correlation energies 
too. Here we can calculate the correlation energy via 
the "ACFD" functional (see e.g. Ref. 21) involving the 
orbital-dependent linear response function xo, and "xc 
kernel" /xc- By way of example, the "PGG"^'' kernel 
directly uses the pair-density to approximate 



f-Kcacr' (^1 ^ ) 



1 



n2Hx 



nc{r)na'(r') 



- 1 



(22) 



It thus captures the ensemble physics at both the LEXX 
and correlation levels via n2Hx- 



V. RESULTS 

In Figure 1 we show correlation- free energies for H, 
Li and Na-like fractional ions calculated in the opti- 
mised EXX and LEXX schemes under the Krieger, Li 
and lafrate^^ (KLI) approximation to the potential in a 
real space code for spherically symmetric systems. Re- 
sults are presented for and ranging from zero to 
one such that / ranges from zero (e.g. Na"'') to two (e.g. 
Na~). The true ensemble EXX energy ijEEXX ^^^j^^g ^j^g 
same, piecewise linear form as (5) but with groundstate 
energies En of the ensemble members (for integer N) re- 
placed by EXX energies Ef^^ from the optimal Slater 
determinant. The sides of the surface plots show the case 
where one electron is integer and the other fractional (or 
integer at the corners) and it is clear that the results for 
the optimised EXX and LEXX are identical as expected. 
In the interior, however, a different picture emerges, with 
the required derivative discontinuities sA f"^ + — 1 be- 
ing absent in the EXX but clearly present in the LEXX. 
The LEXX also varies minimally with / = + /"'' fixed 
(along diagonals perpendicular to the projection), unlike 
the EXX. The slight remaining non-linearity must be ex- 
plained via the implicit dependence of the orbitals on / 



^LEXx^ 1^ i) (Ha) of H, Li and Na ions with fractional oc- 
cupations under EXX (left) and LEXX (right). 

e' - (Ha) EXX e' - e" (Ha) LEXX 




as the energy formula is explicitly linear in /. We are un- 
sure if this is a result of the optimised effective potential 
approach itself, or the KLI approximation thereto. 

The LEXX clearly offers dramatic improvements over 
the EXX in energy calculations. For Li and Na it also 
makes a good approximation to the true EEXX energy 
without resorting to correlation physics. Here the maxi- 
mum variation from EEXX is at most 6mHa for Li and 
Na, significantly smaller than the correlation energies of 
45mHa and 396mHa respectively^^ for the neutral atoms. 
Only for H, where the orbitals of H and H~ differ signif- 
icantly through space, is the difference significant, grow- 
ing to almost 20mHa for / w 1.5, comparable to the H~ 
correlation energy of 42mHa. 

The LEXX was previously used^'^ to generate ground- 
states for correlation energy calculations. We are thus 
able to compare the correlation-free LEXX results from 
that work with benchmark HF energies calculated by 
Chakravorty et alP'^ . To test the validity of the p shell 
LEXX expression (using equations (15)-(17)) we com- 
pared the energies of the first and second row open p 
shell atoms B-F and Al-Cl as these have integer electron 
numbers, but fractional f . For these atoms the LEXX 
energy has a maximum error of < 1.5mHa (for O) and a 
mean average error of just 0.6mHa. To numerical accu- 
racy in our calculations this is close to exact agreement, 
and justifies both the LEXX itself and the KLI approxi- 
mation to the OEP, at least for integer electron number. 

In Figure 2 we show the energy of carbon and fluo- 
rine ions with flve/eight to seven/ten electrons. For il- 
lustrative purposes we show results with (RPA, PGG, 
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FIG. 2. Groundstate energy E{N) of C and F ions under 
the LEXX approach with and without correlation energy in- 
cluded. Electrons are split equally between up and down spin 
N^^N^ = N/2. 



correct treatment of Hx in most cases, with a nonzero 
but very small extra contribution from correlation in the 
PGG case. We aim to further investigate correlation en- 
ergies at fractional occupation in future work. 





exact) and without (LEXX, EEXX) correlation energies 
evaluated in the "ACFD" functional (see e.g. Ref. 21). 
The LEXX is used for the kinetic, external, Hartree and 
exchange energies in all calculations bar EEXX and ex- 
act. Correlation energies are evaluated using the random- 
phase approximation (RPA) and PGG kernel (see Ref. 23 
for technical details). The exact groundstate energy of 
fractional ions of C is given by the piecewise linear func- 
tion E{N) = E§ - (N - 6)/^ for 5 < TV < 6 and 
E{N) = E§-{N- 6)A'^ for 6 < iV < 7 where E§ is the 
groundstate energy of carbon, is its ionisation poten- 
tial and A'~^ its electron affinity (with similar expression 
for F). Energies and ionisation potentials are taken from 
Ref. 28 and affinities from Refs. 29 and 30. The EEXX 
energy is defined in the same way but with Eq, I and A 
replaced by correlation-free EXX values. 

Clearly the LEXX without correlation approximates 
the piecewise linear form, albeit incorrectly predicting 
negative fractional affinities for < 6.75 for C and 
N < 9.60 for F. Including correlation improves things, 
although even here there is a small range with negative 
affinities, at TV < 6.25 for C with the RPA and PGG 
kernels, and N < 9.25 for F with the PGG kernel. It 
is clear that the "LEXX-PGG" (PGG evaluated with 
an LEXX pair-density) is a fairly good approximation 
to the groundstate ensemble energy at all fractions in 
both cases, especially for the positive ions. The deriva- 
tive discontinuity shown here comes entirely from our 



VI. CONCLUSIONS AND FURTHER WORK 

While the discussion here has focused on Fermionic 
systems with non-degenerate frontier orbitals and en- 
sembles constructed around varying electron number, the 
general approach holds true for any non-interacting en- 
semble system. For example in Bosonic systems, or- 
bital SI is not cancelled by exchange terms even for in- 
teger occupation, a situation which favours the present 
type of analysis of the "Hartree" and "exchange" terms. 
Other interesting cases include finite distance dissocia- 
tion, where quantum superpositions of determinants are 
required as well as classical ensembles; and thermal en- 
sembles. 

LEXX physics is also useful beyond the OEP LEXX 
method discussed here. It should be possible to construct 
local density functionals (like the LSDA) from pseudo- 
densities based on the modified exchange and/or Hartree 
pair-density via an approach like that of Ref. 31 or Ref. 
32. This perhaps provides some further justification for 
the success of recent work by Johnson and Contreras- 
Garci'a^^. The LEXX may also have potential uses in 
0(A^)-scaling DFT approaches (see Ref. 33 for a recent 
review) . 

By constructing a density matrix with similar prop- 
erties to the exact ensemble, we were able to develop 
an LEXX formalism yielding an orbital-dependent total 
energy (8) via a pair-density, piecewise linear in the oc- 
cupation factors, and involving ensemble averages of the 
one {9f)£ and two (0f 0J )£ orbital pair factors [see (9)- 
(11)]. This is exemplified for doubly fractional s shells 
in (12)-(13) with ghost-interactions supressed by the cor- 
rection term (14) and with similar expressions for p shells 
discussed in (15)-(16) with additional like-spin correction 
term (17). Using these energy expressions in the OEP 
LEXX functional proposed here gives clearly improved 
resuhs (with £;EEXX < ^lexx < £;EXX^ ^^^^ compared 
with the more common form of EXX, without resorting 
to correlation physics. This suggests that the very notion 
of electron correlation is imprecisely defined for OEP or 
KS systems with fractional occupancy. Using the prop- 
erties of ensembles to create better trial wavefunctions 
and density matrices can be an excellent means of re- 
ducing the workload of the correlation functional in such 
systems. 
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